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Abstract 

Auxiliary bearings are a critical feature of any 
magnetic bearing system. They protect the soft 
iron core of the magnetic bearing during an over- 
load or failure. An auxiliary bearing typically con- 
sists of a rolling element bearing or bushing with a 
clearance gap between the rotor and the inner race 
of the support. The dynamics of such systems can 
be quite complex. It is desired to develop a rotor- 
dynamic model and assess the dynamic behavior 
of a magnetic bearing rotor system which includes 
the effects of auxiliary bearings. Of particular in- 
terest is the effects of introducing sideloading into 
such a system during failure of the magnetic bear- 
ing. A model is developed from an experimental 
test facility and a number of simulation studies 
are performed. These results are presented and 
discussed. 

Nomenclature 

C = damping, N-s/m 

Cl = clearance of auxiliary bearing, m 
D = nominal gap thickness, m 

= external force vector acting on the rotor in 
X direction, N 

Fy = external force vector acting on the rotor in 
Y direction, N 

G = gravitational acceleration, m/s 2 
I = rotor inertia matrix 
i = current, amp 
K = stiffness, N/m 
k = gain value 
L = equivalent circuit length 
M — mass, kg 

N = total number of modes considered 

Nbi = node number at rotor left end 

Nb2 = node number at AMB 

N S 3 = node number at auxiliary bearing 

Nimb = node number at imbalance location 

Nj — number of turns of wire per coil 

Q x = rotor modal coordinate vector in X direction 

Q Y = rotor modal coordinate vector in Y direction 


R H = radius of rotor journal, m. 

R s = radius of auxiliary bearing bore, m. 
t = time, s 
tj = fail time, s 
v = voltage, volt 

v r = relative velocity at auxiliary bearing 
rotor /stator contact point 
Xfl = rotor physical coordinate vector in 
X direction 

Xs = stator physical coordinate vector in 
X direction 

Yk = rotor physical coordinate vector in 

Y direction 

Y s — stator physical coordinate vector in 

Y direction 

a = acceleration, rad /sec 2 
6 = dynamic clearance 

r = 

& = rotor free-free modal displacement matrix 

^ = rotor free-free modal rotation matrix 

if) = imbalance vector 

fi 0 = permeability of free space 

\l t = permeability of silicon steel 

ft = rotor operating speed, rad/s 

r = time constant 

0 = angular position of the shaft 
£ = modal damping coefficient 

Subscripts 
a — first 

B = auxiliary bearing 
6 = second 
C = contact 
g = gain constant 
h = horizontal 

1 = left bearing 

r = rate constant 
v = vertical 
x = x — direction 
y = y — direction 
bi = bias 

1 = top magnet 

2 = right magnet 


1 


3 = bottom magnet 

4 = left magnet 


Superscripts 

a = current amplifier 
c = controller 
l — lead 

p = sensor amplifier 

Introduction 

In recent years, the use of active magnetic bear- 
ings (AMB) for turbomachinery support has been 
an area of interest in both academia and in indus- 
try. Magnetic bearings provide the potential for 
significant improvements over other types of ro- 
tor supports, including elimination of wear, bear- 
ing friction-related energy losses and a means of 
actively suppressing rotor vibration. A critical 
feature of any magnetic bearing supported rotor 
system is auxiliary bearings to protect the soft 
iron core from rotor contact (and subsequent dam- 
age) during an overload or failure of the AMB. If 
the AMB fails partially, i.e. only one coil fails, 
the question arises to which is more advantageous: 
cutting off the AMB entirely or utilizing the re- 
maining coils to aid in predicting the system re- 
sponse. 

The present work is concerned with developing 
an understanding of the response characteristics 
of a rotor system supported by auxiliary bearings 
during failure of the AMB with and without use of 
sideloading. A design for an experimental test rig 
is presented and used as the basis for this study. 
A simulation model is developed in which response 
characteristics are studied and discussed. 

Experimental Test Facility 

The model and control system development 
presented in this study is for an experimental 
test facility. The apparatus consists of a radial 
magnetic bearing that is supporting the right end 
of a rotor. The left end of the rotor is supported 
by a ball bearing suspended in a frame by four 
springs. A photograph and the corresponding 
schematic diagram are shown in Figure 1. 

The apparatus has two basic components: the 
rotor shaft and the magnetic bearing. The shaft is 
made of steel and is 0.0098 meters in diameter and 
0.4572 meters in length. A steel disk of diameter 
7.62 cm, thickness 2.54 cm, and mass 0.810 kg 
is placed at the midpoint of the bearing span. 
Threaded holes on the disk allow for imbalance 
to be added to the system. 



Figure l.a Experimental Test Facility 



Figure l.b Schematic Diagram of Experimental 
Test Facility 


The rotor is driven by a variable speed motor. 
The speed of the rotor can be adjusted by turning 
the appropriate knob on the front panel of the 
control box. Shaft vibration is measured using 
eddy current proximity displacement sensors fixed 
to measure displacement in both the vertical and 
horizontal directions. The displacement signals are 
then sent to a signal analyzer. 

The magnetic bearing is based upon a design 
described in Humphris, et al (1988). The basic 
parameters for the bearing are shown in Table 
1. The bearing consists of four electromagnets 
equally spaced around a soft-iron core. A pho- 
tograph and schematic diagram of the bearing are 
shown in Figure 2. 
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Q y -ftrQ* + c^Q v + * T F y = 0, (1.6) 



Figure 2. a Magnetic Bearing Assembly 



The auxiliary bearing consists of a ball bearing 
with a clearance between the rotor and the stator. 
The governing equations are 


Mb Xb+C bxXb + KbxXb — F Xt0LUX , (l.c) 
MbYb^C B y^B “I" AgyYs = Fy t0UX> (1-d) 


where 


F x,/*n(^61 ) 
Fy t U n (Nbl) 

F x,tmb{Nimb) 

F 

F x,aux(^b3 ) 
F y ( aux( ) 
F grav 

F* 


&N bl (KlxQx + Cix Q x ) 
®N bl (KlyQy + ^yQv) 
Ft 2/ il)§cos6 + a^6 stnO, 
ffifasinO — at/?6 cos6, 
*7V fc3 (FA rcos/3 - F R sin(i) 
^N b3 (F N sinf3 + F rcos(3) 
& t MG 
d>K c (6 - D). 


If tv = 0, then F^ = required static friction 
force. If v r ^ 0 or if required static friction force 
> fi 3 Fat, then F R — fik F/v* 


F x — F Xt /in "I" Fx,»m6 H" F Xt attx F Xl am6 


Figure 2.b Schematic Diagram of Magnetic 
Bearing 


where 


Simulation Model 

A simulation model was developed for the ro- 
tor system described above. The model has three 
principal components — the rotor, magnetic bear- 
ing, and electronics. The governing equations for 
each are shown below. 

The rotor is modelled using the free— free bend- 
ing mode shapes and natural frequencies obtained 
through finite element analysis. The finite element 
code uses 11 stations and the first four modes (two 
rigid body and two flexible modes) are included 
in the simulation model. The rotor equations of 
motion can be expressed in terms of modal coor- 
dinates as 


Sm(/?) \/(Xr - Xs) 2 + (Yh - Y s ) 2 , 

x fl -x 5 

cos(/3) ^/(Xr — Xs ) 2 + (Yr — Y s ) 2 , 

* = ^($Q x -X fr ) 2 + ($Q y -y 6 ) 2 , 
<f> = 1 if S > A 
0 otherwise 
Qx = * T X r , 

Q„ = * T Y r , 


Q x + firQ y + u> 2 Q r + * T F r = 0, ( 1 . 0 ) 
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with , T 

X r = {-Xri,Xr2, • • > 

Y r = {y r l,y r2 ,-..^rm} T . 

(m = total number of nodes) 

The physical displacements of the rotor at the 
bearing and imbalance locations can be obtained 
using the following coordinate transformation: 

N 

X r k — $kiQxi) 

i=1 (k = N bu N h2 ,Nto,Nimb) 

N 

Y r k = ^ &ki Qyi) 

t=l 

The position of the shaft at the magnetic bear- 
ing location is measured using horizontal and verti- 
cal proximity sensors. They are modelled as linear 
first order systems. 


dyl _ k’y - < 
dt t p 


(2 -a) 


Four current amplifiers (one for each coil) sup- 
ply current to the magnetic bearing. They are 
modelled as 


djl 

dt 


*; + — 

T a L T a 1 v 


(5.a) 



dil 

dt 


dll 

dt 


i 2 a + — 

j-a * T a * v 


-~ a * + 


-^r*4 + 


(5.6) 


(5.c) 


(5.d) 


The forces supplied by the magnetic bearing 
to the rotor are assumed to be decoupled in the 
horizontal and vertical directions. The effects of 
flux saturation are included in the model based 
upon the approach described by Lewis (1993). The 
equations for these forces are 


dv p _ k p x - 
dt r p 


( 2 . 6 ) 


The active control system is based upon a 
proportional-derivative (PD) control law and im- 
plemented using analog circuitry. The controller 
is modelled as 


*£ = + Ik, + (3-o) 


^ = w + (*. + (3-*) 


dt 


The control circuitry also includes a lead net- 
work. The governing equations are 


(4 ' a) 


dv 1 


1 , , dv c h 1 j 


— n - ..c i 


(4.6) 


y*,am6 — J 

(2 (D - x Ng3 ) + — ) 2 

Mr 

2± __) 

(2 (D + x Nb3 ) -I- — ) 2 

Mr 

(6.a) 


Fy.amt = Moj4-^t( “ J 

(2(D - 1/ATjjj) + — ) 2 

Mr 

(3 j 

(2(D + j/jvb») + — ) 2 

Mr 

( 6 . 6 ) 


Discussion of Results 

Using the rotordynamic model described in th( 
preceding paragraphs, a study of the combined dy- 
namics of a magnetic bearing supported rotor with 
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Figure 3 Block Diagram of Simulation Model 


Parameter 

Value 

Units 

Ku, K ly 

17,510 

N/m 

Clx> Cly 

2.0 

N-s/m 

K c 

87,557 

N/m 

C c 

0 

N-s/m 

K b 1 

17,510 

N/m 


2.0 

N-s/m 

Mb 

0.033 

kg 

V>6 

3.0 x 10" e 

kg-m 

k p 

7.8 x 10 3 

- 

k 9 

13 

- 

kr 

23 

- 

k a 

0.45 

amp/v 

T p 

1.59 x 10" 4 

- 

^ 

1.8 x 10- a 

- 

rJT\ 

1.59 x 10~ 4 

- 

T' 

5.31 x 10" s 

- 

A 

3.42 x 10" 4 

m 2 

N t 

164 

- 

D 

0.9 x 10“ 3 

m 

Mo 

1.26 x 10 _t! ‘ 

- 

5T“ 

4.2 

volt 



1.9 

volt 


1.0 

volt 


1.9 

volt 


Table 1 Simulation Model Parameters 
auxiliary bearings was performed. The parame- 


ters used in the simulation studies were identified 
from experimental evaluation of the test rig and 
are shown in Table 1. All of the responses illus- 
trated in the figures are for this base parametric 
configuration, unlesss otherwise indicated. The ro- 
tor speed of 200 rad/sec was selected as a reason- 
able value above the first critical speed. 

Figure 4 shows the rotor response without AMB 
failure. The horizontal and vertical responses are 
of the same amplitude and frequency at steady 
state. The response consists of a transient re- 
gion of about 0.5 seconds in length followed by 
steady synchronous oscillations. As expected, the 
responses are similar to those for a conventionally 
supported rotor. 

This study investigates the use of the three re- 
maining coils of the AMB when one fails. By com- 
manding the voltage of the controller, the attrac- 
tive force, or sideloading of the remaining coils can 
be adjusted to help predict where the shaft will hit 
the auxiliary bearing and how it will behave once 
there is contact. This sideloading can be useful 
in providing stability to the system during partial 
failure of the AMB. The test configuration con- 
sists of allowing coil 2 to fail and investigating the 
effects of sideloading on the system. Coils 1 and 
3 are kept in the same electronic configuration as 
if no failure had occurred. However, the voltage 
in coil 4 is set to a value that will pull the shaft 
horizontally over to the auxiliary bearing. Figures 
5 through 7 show the rotor responses for failure 
of a single coil. Of interest in these figures is the 
influence of auxiliary bearing clearance and spin- 
down rate on the response characteristics. For the 
cases examined, the auxiliary bearing is quite stiff 
and the damping levels are quite low. For Fig- 
ures 5. a - 5.c, the right coil of the AMB fails and 
all of the other coils remain active. Since the cou- 
pling between the horizontal and vertical motion is 
primarily through the gyroscopic terms, at steady 
state the vertical responses are little effected by 
the failure of the horizontal coil. It is observed 
that the responses may be quite high as the ro- 
tor bangs around on the auxiliary bearing, even 
with sideloading from the remaining coil. How- 
ever, the sideloading is shown to be quite effective 
in encouraging rotor contact that will dissipate the 
vibrational energy. The response amplitudes tend 
to decrease as the auxiliary bearing clearance is 
decreased. However, a definite trade-off was en- 
countered. If the clearance is too small, the rotor 
may strike the auxiliary bearing during steady- 
state operation if the imbalance is sufficiently high. 
This tends to unnecessarily perturb the rotor and 
shortens the life of the auxiliary bearing. 

Figures 6.a - 6.c and 7.a - 7.c show the rotoi 
behavior for spin-down at various rates. Again 
the response amplitudes tended to be lowest foi 
smaller clearance values. However, an interesting 
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Figure 4 Rotor Response Without AMB Failure 
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Figure 5.b Rotor Response With AMB Failure 
Q = 200rad/sec,CL = ^ 
tf = 0.5, a = 0 
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Figure 5. a Rotor Response With AMB Failure 

° on 

Q = 200rad/sec, C*, = 

tf — 0.5, a = 0 

behavior occurs as the rotor is spun down through 
the critical speeds. In each case, the highest 
amplitude responses are observed for the case 

with a clearance of For the higher and 

lower clearance values, the maximum amplitudes 

on 

are lower. The effective stiffness for the -j- 
case is smaller, so the frequency of the transient 
vibrations is smaller. 

An interesting phenomena is observed for the 
sideloading. The magnetic bearing forces are non- 


Figure 5.c Rotor Response With AMB Failure 
ft = 200rad/sec, Cl — y 

tf = 0.5, a = 0 

linear functions of the dynamic clearance and the 
current. If one coil fails and the remaining coils are 
used to provide sideloading, the natural approach 
might seem to issue a constant command voltage, 
with the objective of having a constant sideload- 
ing force at steady state. However, the inherently 
unstable behavior of the magnetic bearing without 
active control coupled with the vibrational effects 
between the auxiliary bearing and the rotor, pro- 
duces a severely unstable response. A remedy is to 
command a voltage that is proportional to the hor- 
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Position at AMB _ Position at AMB 



Figure 0.a Rotor Response With AMB Failure 
0 on 

ft = 200rad/sec, Cl = -g- 

tf = 0.5, a = —50 



Figure 6.c Rotor Response With AMB Failure 
ft = 200rad/sec, Cl ~ ^ 
tf = 0.5, a ~ -50 
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Figure 0.b Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = y 
tf = 0.5, a = -50 


Figure 7.a Rotor Response With AMB Failure 
ft = 200rad/sec, Cl — ^3^ 
tf = 0.5, a = —100 


izontal shaft displacement. Without an ideal am- 
plifier, any value of desired sideloading force can 
be obtained with such a system. However due to 
the delay in the dynamics from the time constant 
of the power amplifier, unstable behavior may also 
result if the commanded voltage is too large. Fig- 
ure 9 shows the maximum steady-state sideloading 
force that can be commanded with stable behavior 
as a function of power amplifier time constant. We 
see that a sufficiently low time constant will allow 


for virtually any desired level of sideloading. 

Conclusion 

A simulation model has been developed for a 
magnetic bearing supported rotor system with 
auxiliary bearings including frictional effects. The 
model has been described in detail with exper- 
imentally obtained model parameters. The re- 
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Figure 7.b Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = y 
tf = 0.5, a — —100 


Figure 8 Rotor Natural Frequencies as 
a Function of Rotor Speed 
(for the parameters of Table 1) 



Figure 7.c Rotor Response With AMB Failure 
ft = 200rad/ sec, Cl = y 
tf = 0.5, a = —100 

sponse characteristics for a variety of auxiliary 
bearing clearances and spin-down rates are pre- 
sented and discussed, including the effects of side- 
loading. In addition, some guidelines are given for 
the selection of appropriate levels of sideloading. 
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Appendix — Linear Current and Position 
Stiffnesses 

For design purposes, a linearized model for the 
magnetic bearing forces is needed. 

The current stiffnesses are: 


„i n 0 Ak\v\' Nt 
Kl ~ 2D* 


ii 0 Ak$v\'Nl 

* 2 ~ 2 D 2 

p<>Ak% v$*Nt 

W 2 


(7.a) 

(7.6) 

(7.c) 


K\ = 


p 0 Ak 


2 U 2 


bi N$ 


2D 2 

The position stiffnesses are: 


K( = 


K? = 


PoAjk^fN]. 

2D 3 


PoAjk^fN^ 

2D 3 


(7 -d) 

(8.o) 

( 8 . 6 ) 
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